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Abstract 



U(l) family symmetries have led to successful predictions of the fermion mass spectrum 
and the mixing angles of the hadronic sector. In the context of the supersymmetric unified 
theories, they further imply a non-trivial mass structure for the scalar partners, giving 
rise to new sources of flavor violation. In the present work, lepton flavor non- conserving 
processes are examined in the context of the minimal supersymmetric standard model 
augmented by a [/(l)-family symmetry. We calculate the mixing effects on the fi — > 
and r — > /i7 rare decays. All supersymmetric scalar masses involved in the processes are 
determined at low energies using two loop renormalization group analysis and threshold 
corrections. Further, various novel effects are considered and found to have important 
impact on the branching ratios. Thus, a rather interesting result is that when the see-saw 
mechanism is applied in the 12 x 12 sneutrino mass matrix, the mixing effects of the 
Dirac matrix in the effective light sneutrino sector are canceled at first order. In this class 
of models and for the case that soft term mixing is already present at the GUT scale, 
r — > /i7 decays are mostly expected to arise at rates significantly smaller than the current 
experimental limits. On the other hand, the /i — > rare decays impose important bounds 
on the model parameters, particularly on the supersymmetric scalar mass spectrum. In 
the absence of soft term mixing at high energies, the predicted branching ratios for rare 
decays expected, well below the experimental bounds. 
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1 Introduction 



In the last years there has been a lot of interest in lepton flavor violation, which can 
be a powerful tool in any attempt to classify different extensions of the Standard Model 
(SM) [1-4]. Indeed, there are various ways to enlarge the known particle spectrum in a 
manner that lepton number violation is generated. The new states may be right-handed 
neutrinos, additional vector-like heavy fermions, supersymmetric particles, or even all the 
above together as this is the usual case in supersymmetric grand unified theories (SUSY 
GUTS) and their string versions. 

A well known result in the context of the non-supersymmetric standard model is the 
conservation of lepton flavor in the case of zero neutrino masses, while in the case of 
massive, non-degenerate light neutrinos, the amount of lepton flavor violation is propor- 
tional to the factor r] u = Am^/M^ ||. However, non-observation of double beta decay 
as well as other well known neutrino data, imply severe bounds on the light neutrino 
mass-squared differences Ami, leading to a high suppression of the ratio rj v and therefore 
to all flavor violating processes depending on it. An analogous smaller suppression occurs 
if the process is mediated by heavy neutrinos. When supersymmetry enters in the game, 
the whole scene changes completely. Even in the absence of right handed neutrinos, flavor 
violations could occur via the exchange of supersymmetric particles. A large number of 
new parameters (sparticle masses, mixing angles, e.t.c.) appear in the calculations, there- 
fore enlarging the number of possible decays, while the predicted branching ratios are now 
comparable to the present experimental bounds. Moreover, the soft breaking terms may 
violate the individual lepton numbers L e , L M and L T by large amounts. Therefore, in the 
context of unification and low energy phenomenology scenarios, flavor- violating processes 
can provide useful constraints on the parameter space of a given model. Even in models 
with universal initial conditions for the scalar masses, the renormalization group runs of 
the slepton masses will give rise to lepton flavor violation that can be significant. 

Thus, there are various sources of Lepton Flavor Violation (LFV) in SUSY models that 
may be tested in future searches ||. Then, it is interesting to investigate them in the 
framework of flavor symmetries, which will give us a guideline on the amount of flavor 
violation that we may expect at the GUT scale. Indeed: 

• 1) Unified supersymmetric theories with U(l) family symmetries generate success- 
fully the observed low energy hierarchy of the fermion mass spectrum and the quark 
mixing in terms of a minimum number of arbitrary parameters at the unification 
scale. In addition, a non-diagonal structure of the sparticle mass matrices comes 
out as a prediction of the theory. Rare processes are sensitive to the scalar mass 
matrix structure and non-diagonality of the latter usually, in a basis in which the 
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fermions are diagonal, leads to hard violation of flavor. 

• 2) These flavor violating effects are enhanced in particular when the higgs vacuum 
expectation value (vev) ratio (tan j3) is large. In fact, many models based on a single 
unified gauge group (like SO (10) theory) predict equality of the top and bottom 
Yukawa couplings at the GUT scale, and therefore a large value of tan j3 is implied. 
As a result, the 6x6 structure of the slepton mass matrix enhances further the 
lepton mixing effects. 

• 3) When right-handed neutrinos enter in the model, the theory faces another chal- 
lenge. Firstly, Dirac mass matrices arise of the order of the up-quark masses. In 
practice the majority of unified theories imply a Dirac matrix at the unification scale 
equal to the up quark matrix. Charged leptons and neutrinos are no longer diagonal 
in the same basis and a leptonic mixing matrix, similar to the Cabbibo Kobayashi 
Maskawa matrix Vckm for the quarks, is unavoidable. Secondly, the Dirac mass 
matrix itself has an even more intriguing role, since, due to renormalization effects 
on the supersymmetric scalar sector, it modifies the slepton mass matrix at low en- 
ergies. Moreover, it enters in the construction of the 12 x 12- sneutrino mass matrix 
which in principle would have the potential to give rise to additional flavor- violating 
effects. Nevertheless, since only the effective light sneutrino mass matrix is relevant 
in the calculation, we will show in this work that the mo effects are canceled at first 
order, when the see-saw mechanism is applied. 

• 4) The contributions of the trilinear scalar mass parameter of the potential, i.e, 
A-term contributions, are also discussed in this work in some detail. Due to the 
absence of direct experimental information about its values, sometimes its effects 
are not considered at all. However, imposing even a zero value for initial condition 
at the GUT scale, its low energy value will not be unmarked. Renormalization group 
running effects will drive its value to magnitudes comparable with those of scalar 
masses, whereas its contribution to the branching ratios is of particular importance. 

In this work, we calculate the branching ratios for lepton flavor violating decays in 
models with abelian flavor symmetries which emulate the situation in many string con- 
structions. We concentrate in particular in the class of models for fermion masses with 
one family symmetry which is the simplest possibility and have been firstly proposed in 
0. This calculation offers an important test with regard to the viability of these particu- 
lar models, but also gives further insight on the general predictions of family symmetries. 
The motivation for their introduction is to explain the observed fermion mass hierarchy 
yet they have further consequences as is the case of flavor violations. The lepton sector is 
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ideal for such checks (in the quark sector, large uncertainties may enter in the calculations, 
due to poor knowledge of hadronic matrix elements). 

The analysis follows the lines of ||, where some first estimates for flavor- violation in 
this class of models have been presented. This paper is organized as follows: 

In section two we briefly analyze the basic features of a generalized class of SUSY 
models with one C/(l)-family symmetry. We give the forms of fermion and scalar mass 
matrices and discuss their role in flavor violating parameters. 

In section three, we derive the 12 x 12 sneutrino mass matrix and show that the effective 
3x3 light sneutrino matrix entering the flavor violating decays does not involve Dirac 
neutrino mass contributions. 

In sections four and five we give the loop calculations for the various amplitudes and 
analyze the procedure for our numerical investigations. 

In section six we present the results treating separately the cases with and without 
scalar mass mixings at Mqut- Finally, in section 7 we present our conclusions. 

2 Mass matrices 

As has been stressed in the introduction, one of the main advantages of the [/(l)-family 
symmetries (U(l) /) is the determination of the hierarchy of the mass spectrum and mix- 
ing angles and the prediction of many other parameters of the minimal supersymmetric 
standard model (including the scalar matrices) with only a minimal set of arbitrary pa- 
rameters. In the simplest case of only one U(l)f symmetry, the fermion mass hierarchy 
is successfully obtained using an additional singlet field which develops a vev one order 
of magnitude below the string scale Mjj. Throughout our calculations, we will further 
assume the existence of a grand unified symmetry at a scale M GUT without specifying 
the gauge group. (We have in mind models with intermediate gauge symmetry groups|| 
where the gauge couplings run together from M\j down to Mqut-) Below the unification 
scale, only the minimal supersymmetric spectrum is assumed, therefore the unification 
point will be taken at ~ 10 16 GeV. In this scheme, the low energy parameters involved in 
our subsequent calculations depend on M GUT , the common value of the gauge coupling 
gcur and the ratio of the singlet vev over the string scale, e ~ ((f)) /My. When the U(l)f 
symmetry is exact, only the third generation has a Yukawa term in the superpotential, 
whilst all mixing angles are zero and lighter generations remain massless and uncoupled. 
When <j) acquires a vev, the [/(l)j-symmetry is broken and mass terms fill in the rest of 
the mass matrix entries with Yukawa terms suppressed by powers of the ratio (j)/Mjj. In 
simple supergravity models, a similar situation occurs also in the scalar sector; at tree 
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level, however, there appear three instead of one, diagonal mass terms, one for each gen- 
eration. Mixing (off-diagonal) terms in the scalar sector shows up when higher NR-terms 
are included in the Kahler potential. Thus, U(l) /-symmetries imply also a non-trivial 
structure for the corresponding scalar mass matrices. This additional structure may be 
responsible for new hard flavor violations. The existing experimental bounds on flavor- 
violating processes are therefore going to give us an indication as to which structures are 
viable for the scalar mass matrices. 

Which are the elements that determine the form of the scalar mass matrices? Clearly, 
the charges of the various fields under the flavor symmetry will play a dominant role. 
Moreover, as we discussed in the previous paragraph, the mass and mixing hierarchies 
depend on expansion parameters that are generated when singlet fields acquire vev's. 
What these vev's can be, depends on the flat directions of a given theory. Once the flat 
directions and the U(l)f charges of a particular model have been fixed, the scalar mass 
matrix structure may be easily computed through the Kahler function Q = JC + log | W| 2 
where W is the superpotential and /C has the general form 

K = - \og{S + S*)-J2 h n log(T n + T* n ) + Z tj . (T„, T^Q.Q* + ■■■ (1) 

with Qi being the matter fields, S the dilaton, whereas T n are the other moduli fields. 
The scalar mass matrices are determined by Z^* and W. The form of the Z^* function is 
dictated by the modular symmetries and depends on the moduli and the modular weights 
of the fields. Thus, at the tree level, the diagonal terms are the only non-zero entries in 
the scalar mass matrices. Higher order terms allowed by the symmetries of the specific 
model fill in the non-diagonal entries. 

The lepton Yukawa interactions which appear in the superpotential in the presence of 
the right handed neutrino are 

We P = e^XJH, + N C \ D £H 2 + \ nX N c N c . (2) 

Here t is the left lepton doublet, e c is the right singlet charged lepton, N c is the right- 
handed (RH) neutrino and X e ,D,N represent Yukawa coupling matrices in flavor space. 
Also, Hi and H 2 are higgs doublets and x stands for an effective singlet which may 
acquire a vev at a large scale. 

In addition, soft supersymmetry breaking terms generate mass matrices for the charged 
slepton fields, denoted by m,£, rh eR - Denoting the various fields collectively with z, in 
supergravity the scalar potential is given by 

V = e G{z) (dGjj-Gj - 3) + I D | 2 (3) 
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where | D \ 2 represents the contribution of the D-terms in the potential. Also, with Gj, 
we denote the derivatives of G with respect to the fields z T , i.e., 

Gj = ±V,W (4) 

where T>iW = djW + W<9//C is the Kahler derivative. Writing explicitly the various fields, 
in the low energy limit one encounters the following scalar mass terms in V 

V = m\H{H x + m 2 2 H* 2 H 2 + m?g*g + m\ c u c *u c + m%d c *d c 

+ mjti + m%~e c *e c + m% c N c *N c 

+ {e ab (mlH«H b 2 + A u \ u q a u c H b 2 + A d \ d q a d c Hl 

+ A^t-^Hl + A v \ D tN c H b ) + A N X NX N C N C }) + h.c. + q.t. + • • • , (5) 

In the above equation a and 6 are SU(2) indices. The dots stand for scalar trilinear 
F-terms and q.t. denotes quartic terms in the scalar fields. Also terms proportional to B 
parameter are not shown explicitly here. The Higgs mass terms contain two contributions: 
the first arising from the superpotential and the second from the soft supersymmetry 
breaking parameters: mf = fi 2 + m 2 ?., with i = 1,2. 

We present here the relevant mass matrices of a model whose successful fermion mass 
hierarchy is predicted by U(l)f symmetries. We work in the low energy effective model 
based on the 577(3) x SU(2) x U{1) gauge group with an additional U(l)f symmetry 0. 
After the implementation of this symmetry, the fermion matrix for charged leptons in this 
model is given by 

/ ^2\a+b\ ~\a\ g|a+6| \ 
me ~ g|a| g2|6| mr (6) 

v ~ e W+b\ ~ e \b\ x ^ 

where the parameter e is some power of the singlet vev scaled by the unification mass, 
while a, b are certain combinations of the lepton and quark U(l) /-charges. Order one 
parameters cy in front of the various entries (not calculable in this simple model) are 
assumed, to reproduce the fermion mass relations after renormalization group running. 
These parameters c^- are usually left unspecified, here however their exact values are 
necessary for a reliable calculation of the lepton violating processes. 

A phenomenologically viable choice for the charges is to take a = 3 and 6=1. A 
successful lepton mass hierarchy in this case is obtained for the choice e = 0.23. In this 
case, a possible choice of the coefficients is given by cu = C21 = 0.4, C22 = 2.2, with the 
rest of the coefficients being unity. We point out that the choice of the parameter b is not 
completely determined by the lepton mass texture. In fact there is a second possibility, 
with 6=1/2 where the fermion mass hierarchy is also consistent with the low energy data. 
On the contrary, the choice of 6, as well as the choice of coefficients, will have a significant 
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impact on the magnitude of the rare processes. Thus, flavor violation is a powerful 'tool' 
and an invaluable criterion of the viability of a certain choice. A detailed discussion on 
this question is one of the main points of the present analysis and will be given after the 
calculations on the branching ratios will be presented in the subsequent sections. 

The Dirac mass matrix in the above model has a similar structure. Due to the simple 
U(l) structure of the theory, the powers appearing in its entries are the same as the lepton 
mass matrix; however, the expansion parameter is in general different [[[(J . Thus, its form 
is given by 

/ £ 2|a+b| e |o| £ |a+b| \ 

m^w eW m top (7) 

v e |a+6| e |6| 1 j 

The choice of charges a = 3, b = 1 allows to identify the Dirac mass matrix with the up- 
quark mass matrix. Again, order one coefficients have to be introduced in the quark mass 
matrix in order to obtain consistency with the experimentally determined masses. We 
denote here the corresponding coefficients multiplying the entries of (0) with dij. A choice 
of coefficients leading to correct up-quark masses is obtained for dyi = di\ = -5, d^2 = 
d 2 s = 1.5, with the rest of the coefficients being unity. In the case of the up-quark matrix 
a second expansion parameter]?! is introduced with the value e = .053. 

The RH-Majorana mass matrix is constructed from terms of the form xiV c iV c where 
X is an effective singlet. In GUT models, \ is a combination of scalar (Higgs) fields. 
Obviously, the structure of Mm depends on the origin of the singlet x as we U as its 
charge Hl0|, |XTf| . Thus, in a class of models this singlet may arise from the combination 
N C N C where N c is the scalar component of the RH-antineutrino supermultiplet. There 
are therefore various structures of the Majorana matrix, depending on the specific choice 
of the x-charge. However, as we are going to show analytically, due to cancellations, the 
results are not sensitive to the structure of Mn Q. 

The scalar mass matrices of this model are built using the potential mentioned in the 
beginning of this section. In particular, for the sleptons we obtain at the GUT scale 

/ 1 gl«+26| l\a+b\ \ 

~\a+2b\ i ~\b\ 

~\a+b\ ~\b\ 1 



m 3/2 



At this point, we have determined all the necessary ingredients in order to build the 

basic quantity which determines the flavor violations in SUSY theories. This is the 6x6 

lr This, we checked numerically, by picking particular forms of Afjv, which also fit the light neutrino 
data. For example, a zero singlet charge leads to a M/v form similar to that of the Dirac mass matrix. A 
singlet charge Q = — 1 leads to an interesting form with large mixing in the 2-3 generations and a two- fold 
degeneracy, suggesting a solution of the atmospheric neutrino puzzle through the v ^ v T oscillations Q. 
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slepton mass matrix which takes the form 

rhj ((Ae + fi tan. (3)me) 

(Ae + fi tan (3)m e m\ R 

mj er are scale dependent, their value being given by the renormalization group equations 
Their initial conditions at Mqut are defined by the matrices (|8|) and have a significant 
effect in flavor violations. In addition, as it has been extensively discussed in the litera- 
ture, due to RGE running slepton mass matrices receive corrections proportional to the 
Dirac mass matrix. Charged lepton and Dirac mass matrices are not simultaneously di- 
agonalized; thus in the basis where mi is diagonal, the 3x3 slepton mass matrix acquires 
a non-diagonal contribution of the form 

6 ^ ^ J^^ 3 + a ^ ln ~^^ X D X D m l/2 (10) 

where \d is the Dirac Yukawa coupling and the proportionality factor depends on the 
scalar mass parameters squared while the parameter a is related to the trilinear mass 
parameter A\ = am 3 / 2 . Thus, even if one starts with a diagonal slepton mass matrix, 
there are important off-diagonal contributions due to the existence of N c . This effect, 
however, will prove less important in theories where the scalar mass matrix textures are 
also determined by U(l) symmetries. 




3 Sneutrino mass matrix 

The sneutrino mass matrix is also determined similarly. It is a 12 x 12 structure given in 
terms of the 3x3 Dirac, Majorana and slepton mass matrices. It is generally expected 
that -as in the case of charged sleptons- the Dirac term induces considerable mixing 
effects. We will show here that this is not the case in the sneutrino mass matrix. 

This 12 x 12 matrix is rather complicated and not easy to handle. Vastly different scales 
are involved and numerical investigations should be carried out with great care. Its form 
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is as follows: 



v 



N c 



m* D M T (A* u + ficot/3) m* D 



M*m T D 



(A u + jicot[3) mo 



(A* u + ficot(3) m 2 N + M*M 7 



mj) (A u + ficotfi) 



m D M+ 



A N *M* 



A N M m 2 N + MM + 

+m D m~] D 



One can construct an effective 6x6 matrix for the light sector, by applying matrix 
perturbation theory, similar to the see-saw mechanism. The result up to second order is: 



[m r , 



eff 



(m D M- 



{(2A U + A N ) + 2/icot/3)- 
(m D M~ 



_1 ™d) 



V 



{{2A U + A N ) +2fjicotp)- 
(moM- 1 ™^)* 



m 



(A v + ncot/3)(A v -2A N )- 
(m D M~ 2 m) D ) 



J 



;i2) 

The first and second order terms are obtained assuming all parameters as real and the 
A— matrices proportional to the identity. Notice that the second order terms along the 
diagonal can be neglected. The first order off-diagonal terms must be retained, since they 
lead to complete mixing of the pairwise degenerate states. This, however, does not affect 
the flavor-violating branching ratio. 

The simplicity of this result is rather astonishing. We note that after the 'see-saw' 
mechanism is applied, the Dirac neutrino mass matrix contribution in the effective light 
snetrino mass sector is essensially negligible. Moreover, there is an additional benefit, 
since the complication of the initial 12 x 12 mass matrix can now be avoided. A direct 
numerical calculation of mass eigenstates and mixing angles would be a hard task, due to 
the vastly different scales. 



4 Amplitudes for flavor violating processes 

Figure 1 shows the one-loop diagrams relevant to the /x — > process. The corresponding 
r — > /!7-decay is represented by an analogous set of graphs. There are also box-diagrams 
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contributing to this process; they are however relatively suppressed. 

The electromagnetic current operator between two lepton states Zj and lj is given in 
general by 

Tx = (k{p-q)\JxHp)) 

= Ui{p- q){m 3 ia xp q P (A L M P L + A m Pr) } Uj (p) (13) 

where q is the photon momentum. The Am's have contributions from neutralino-charged 
slepton (n) and chargino-sneutrino (c) exchange 

aL,R — aL,R . aL,R 
M — A M{n) + %(c) 

The amplitude of the process is then proportional to T\€ X where e A is the photon polar- 
ization vector. An easy way to determine the loop momentum integral contribution to 
the Am's is to search, in the corresponding diagram, for terms of the form (p ■ e) and 
make the replacement 2(p ■ e) — > lax^e^. Defining the ratio x = M 2 /m 2 , where M is 
the chargino (neutralino) mass and m the sneutrino (charged slepton) mass, the following 
functions appear in the Am term 

Au(n) '■ q(i- x ) a (1 — 6x + 3a; 2 + 2a; 3 — 6a; 2 log a;) (L — L amplitude) 
(1-1)3 (1 — x 2 + 2x log x)y/x (L — R amplitude) 

Am{c) '■ 6(i- x )<± + 3a; — 6a; 2 + x 3 + 6x log a;) (L — L amplitude) 



(15) 



(l-x) 



^■(—3 + Ax — x 2 — 2\ogx)y/x (L — R amplitude) 



where is the mass of the lj lepton. 

Notice in the L-L amplitudes the lack of terms proportional to the gaugino mass M 
which cancel. The Branching Ratio {BR) of the decay lj — > I4 + 7 is given by 

BR(lj - la) = ^ {{A L M f + (Am)' 



5 Inputs and Procedure 

The branching ratio formulae for the fi — > c-f and r — > /ry decays involve the masses 
of most of the supersymmetric particles. It is important therefore for any given set 
of GUT parameters to know precisely all masses and the other low energy parameters. 
In the present work, this is obtained by numerical integration of the renormalization 
group equations of the MSSM with right handed neutrinos. The renormalization-group 
equations can be found in many papers (see for example [|T3|] ) . 
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We evaluate the coupling constants, using renormalizat ion-group equations at two loops. 
Threshold effects are also taken into account, by decoupling every sparticle at the scale 
of its running mass Q = rrii(Q). Below the scale m t , we use the SM beta functions.. 

Our analysis uses as input values the unified coupling constant (Xq, at the GUT scale, 
the third generation Yukawa couplings At, A&_ r , the common scalar diagonal scalar masses 
777,3/2, the gaugino mass m^, the (effective) Higgs bilinear coupling //, and the ratio of 
the Higgs vev's described by tan/3 and the flavor-symmetric soft-breaking parameter Aq. 

Our integration procedure consists on iterative runs of the renormalizat ion-group equa- 
tions from M GUT to low energies and back, for every set of input parameters TO1/2, 7723/2, A 
and tan/3, until agreement with experimental data is achieved. The values for a G and 
Mgut are obtained consistently with a em , a 3 and sin 2 9w at m>z- Supersymmetric correc- 



tions to sin Ow are also considered [[U|] . 

A similar procedure is followed to obtain the GUT values for the third generation 
Yukawa couplings. In particular, the GUT value for X t is adjusted by requiring the top 
physical mass to be m t = 175±5GEV. Similarly, we obtain the value of the unified A;,_ T , 
by requiring the correct prediction for m T = 1.778GEV. In all the cases that we analyzed, 
values for mt, consistent with experiments are found once QED and QCD corrections are 
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The value of the \x parameter (up to its sign) can be expressed in terms of the other input 
parameters by means of symmetry breaking conditions. To this end, we use the semi- 
analytic formulae including one loop corrections, as given in ref. |L6||. Finally, the chargino 
and neutralino masses are obtained by diagonalization of the 4x4 neutral and the 2x2 
charged matrices as described in |T7[. The RG evolution of the eigenstates, starting from 
universal initial conditions at the GUT scale, is properly taken into account. 

We then explore the values of the BR for all significant values of the input parameters 
7771/2,7773/2, A and tan/3. 

If we consider common scalar masses and trilinear terms at the GUT scale, leptons and 
sleptons will be diagonal in the same superfield basis. However, due to the presence of (a) 
the non diagonal GUT terms A at the GUT scale, and (b) the appearance of Ad in the 
RG equations, the lepton Yukawa matrix and the slepton mass matrix can not be brought 
simultaneously to a diagonal form at the scale of the heavy Majorana masses. Therefore, 
lepton number will be violated by the one loop diagrams of fig. 1. 

We define the unitary matrices diagonalizing the Yukawa mass textures A73 and A e , as 
follows 

X 6 D = Tl\ D T L (16) 
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K = v£\ e v L (17) 

Here, the index 5 indicates a diagonal form. Then, the mixing matrix K in the lepton 
sector, defined in analogy to Vckm is given by the product 

K = T[V L (18) 

The charged slepton masses are obtained by numerical diagonalization of the 6x6 matrix 



>>> 2 I m l L m l R 1 (19) 

m RL m RR 



where all entries are 3x3 matrices in the flavor space. In the superfield basis where A e 
is diagonal, it is convenient for later use to write the 3x3 entries of (|19|) in the form: 

m\ L = (ml) 2 + Sm 2 N + A L + m] + Mf(~ - sin 2 8 w )cos2(3 (20) 



m RR = (m d e ~ R ) 2 + A R + mf - M 2 z sin 2 6 w cos2(3 (21) 

m 2 RL = (A s e + 5A e + ntan/3)mi (22) 

m LR = m RL ( 23 ) 

Each component above has a different origin and gives an independent contribution in 
the Branching Ratios. We further wish to emphasize the following: 



K) 2 , (m s e - R ) 2 , 



A s e denote the scalar diagonal contribution of the corresponding ma- 
trices; their entries are obtained by numerical integration of the RG equations as 
described before. We consider m\j 2 as the common initial condition for the masses 
at the GUT scale, while the trilinear terms scale as am 3 / 2 - Since in the RGEs we 
consider only third generation Yukawa couplings and common initial conditions at 
the GUT scale for the soft masses, our treatment is equivalent to working in super- 
field basis, such that: (i) is diagonal from the GUT scale to the intermediate 
scale and (ii) A e is diagonal from the intermediate scale to low energies. The change 
of bases will produce a shift in the diagonal elements of the soft mass matrices at the 
GUT and at the intermediate scale. This effect is negligible (less than one percent). 

5m 2 N and 5Ai stand for the off-diagonal terms which appear due to the fact that Ad 
and A e may not be diagonalized simultaneously. The intermediate scale that enters 
in the calculation (which is the mass scale for the neutral Majorana field Mjv) is 
defined by demanding that neutrino masses ~ leV are generated via the "see-saw" 
mechanism. This sets the scale to be around the value 10 13 GeV. Then, the 
following values are obtained: 

Sm 2 N = K x mf(m N ) K\ nondiagonal (24) 

5Ai = VlA^TUnWI I nondiagonal (25) 
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• The following values for Ax and A# are defined at the GUT scale: 

A L = VlAV L (26) 
A R = V^AVr (27) 

The effective 3x3 sneutrino mass matrix squared has the same form as the ra\ L part of 
the 6x6 charged slepton one, with the difference that now Dirac masses are absent (in 
consistency with what we have shown in the analysis of the 12 x 12 sneutrino matrix). 
Thus, 

ml = {ml) 2 + 6m 2 N + A L + ^M 2 z cos2(3 (28) 

It is illustrative to write our results as an approximate function of the input parameters. 
Below we give the numerical range for the sneutrino and the A parameter as these are 
defined in 



/ 



5m 2 N « 



(4.2-6.3) x 10" 



(4.2 - 6.3) x 10 ■ u 
v (2.3-3.3) x 10~ 4 (0.7-1.1 



x 10 



-2 



(2.3-3.3) x 10~ 4 \ 
(0.7-1.1) x 10~ 2 

/ 



(3 + a 2 )ml /2 (29) 



5 A, 



( (1.2 - 1.7) x 10 -4 (5.2-7.3) x 10" 4 \ 

(1.2 - 1.7) x 10" 4 (1.0 - 1.4) x 10" 2 

v (5.2 - 7.3) x 10~ 4 (1.0 - 1.4) x lO -2 j 



A 



(30) 



In the last two equations, the ranges in parentheses correspond to tan (3 values between 
14 and 3 (larger contributions are obtained for smaller values of tan /3). 



6 Results 

We have seen that, when [/(l)-family symmetries are taken into account non-diagonality 
in the mass matrices is generic in both the fermion and the scalar sector. This may 
generate unacceptably large flavor-violating effects. It is possible that cyclic permuta- 
tion symmetries between generations and universal anomalous [/(l)-factors may prevent 



mixing effects in the supersymmetric mass matrices ||18|| . In our results, we are consider- 
ing separately two distinct cases: First, we will consider the case where the scalar mass 
matrices are protected from mixing effects by some kind of symmetry not affecting the 
fermion mass sector. Second, we will allow mixing effects in both sectors, and recalculate 
the branching ratios and the new bounds obtained on the sparticle spectrum. 
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1). Case without scalar mass mixing at the GUT scale, A = 0. 

We start with the process \i — *> cy, in the absence of mixings at the GUT scale. Let 
us denote by (n) the contributions from the neutralino-charged slepton exchange, and by 
(c) the ones from the chargino-sneutrino. Then, 

aL,R aL,R i aL,R 

M — A M(n) "+" AfcT(c)- \ 6L ) 

The various amplitudes appear in Fig.2. As we can see, the two contributions to Afj 
(dashed lines) are of the same order of magnitude and opposite signs, while their magni- 
tude decreases with tan/3. Around a certain value of m 3 / 2 , there is a partial cancellation 
of both amplitudes, leading to a decrease of the expected BR(fj, — > e + 7) , as we can see 
in Figs.0, 1 and § 

The contribution to comes almost exclusively from A^^, since the chargino ex- 
change contribution A\ { ^ (Feynman diagram b in Figure 1) arises due to Yukawa in- 
teractions and is about three orders of magnitude smaller than the other contributing 
amplitudes. It is important to make this remark, since when there are no mixings in the 
right-handed slepton masses, AL n i become relevant due to the presence of non-diagonal 
mixings from the trilinear terms 6 A (otherwise, these contributions are of the same order 
as Ajtf, c -)). The effect of 5 A can be seen clearly in Fig. [|. For the initial condition Ao = 0, 
5 A is significantly suppressed, and hence a dramatic decrease of the BR(fi — > e + 7) is 
observed. On the other hand, the cases with A ^ have a remarkable difference with 
the previous one. All curves now are smoother while there are no particular m 3 / 2 -values 
where BR exhibit large suppression. 

Although we start in our case with universal soft masses at the GUT scale, there is an 



analogy with the situation of the SU (5) model discussed elsewhere |L9], |20| in the following 
sense. Assuming we are in a basis where Yukawa matrices are diagonal, the renormal- 
ization of the universal soft mass terms from Mqut down to the RH-neutrino mass scale 
Mtv will split mass parameters of different generations. As a result, flavor universality in 
the scalar sector at is lost. In the SU(5) case, the deviations from universality arise 
due to the renormalization group running from the Planck scale Mpi down to the SU(5) 
scale Mqut- Lepton flavor violation diagrams arise due to the non-universality of the 
right-handed slepton masses at the GUT scale. In contrast to the case we present here, 
the main contributions to the BR{^i — > cy) arise from the neutralino exchange diagram 
amplitudes Aj^^, while there is a similar cancellation to the one described in the above 
paragraph due to the fact that the two main contributions enter in the calculation with 
opposite sign. f\ In 5*0(10) unified models EH], LFV arises due to the non-universality of 



2 In | fL9| , the LL and LR diagrams are discussed separately, while in ours both contributions are 
included in the single neutralino exchange. 
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both right and left sleptons, and hence this cancellation does not take place. 

We may further compare our results directly with the results of |Z2|| , where a similar 
effect is observed, for the string-embedded version ALR of the Pati-Salam model [24]] |. 
In this work, large values for tan/3 are considered and the contribution arising from the 
chargino exchange diagram is the dominant one. Moreover, in this analysis of vl terms, 
Dirac masses in addition to the soft left lepton masses are also included. For a certain 
value of m.3/2, the renormalization-group effects on the soft masses are canceled by the 
Dirac masses, hence the masses of become universal and the contribution from the 
chargino diagram vanishes. Then, the BR(fi — > e + 7) decreases dramatically. Although 
this effect agrees with our results, we would like to emphasize that the small value of the 
BR for certain values of m 3 / 2 , arises due to a cancellation of two contributing amplitudes. 
Note that the contribution of the Dirac matrix to the sneutrino masses is absent in our 
analysis, as we have shown in our treatment of the full neutral scalar mass. 

Figs.|3], ^ and |5|, show the effects of the changes of the input parameters in the total 
BR(^i — > e + 7) . The changes in the non-diagonal elements of the scalar matrices can be 
induced from the approximate formulas that we presented, however the behavior of the 
BR is not correlated in all cases to the increase in the mixing. More precisely: 



• Fig.|3| shows the increase of the BR as the gaugino masses decrease. We can see that 
for a fixed value of tan/3 = 7, and values of m 1( / 2 leading to a SUSY mass spectrum 
inside the experimental limits, the predicted values for the BR are two orders of 
magnitude lower than the experimental bounds. 

• Fig.|] shows the change of the predicted BR for fixed values of mm. Here, we can 
see an increase of the BR with tan/3 (solid lines). The dashed lines are chosen in a 
way that maximal BR's are obtained; still we can see that the BR's stay below the 
experimental bounds. 

• Fig.[5| shows the effect of Aq in the calculation. In this case, the behavior of the BR 
is directly correlated to the increase in the scalar mixings with A . The sign of A Q 
has very little influence in the total result, since it is the square of the A— parameter 
that enters in the relevant renormalization-group equations. 



2). Case with scalar mass mixing at the GUT scale, A^O. 

Let us now proceed to discuss in detail what happens in the case that mixing of soft 
terms (arising through Al,r 7^ 0) is already present at the GUT scale. Here we should 
note that the contribution from A^ 7^ (which is independent of tan/? and m^), 

3 The phenomenology of these string inspired models has been discussed in GS]. 
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is much bigger than 5m 2 N . This implies that there will be a dramatic increase of all 
the Branching Ratios and all previously noted effects have to be reconsidered. Note for 
example that, since m RR is now non-diagonal, additional amplitude contributions destroy 
the cancellations that we observed in the previous case. Fig. 6 shows that the absolute 
value of A\ { is now bigger than v4f f , (the contribution from Aj^(c) to is again three 
orders of magnitude smaller than the others, since it involves Yukawa couplings). 

The modifications in the case that mixing of soft-terms occurs at the GUT scale, are 
presented in Figures [7], |8|. Changes of the BR with A Q are less relevant in this case 
since, as we mentioned before, the dominant mixing terms in the soft mass matrices are 
independent of it. We have used A = — 1.5m 3 / 2 in all the calculations of the rest of the 
section. 

Figure [7| indicates the variation of the BR as a function of m 3 / 2 , for m\/2 = 200, 300, 500 
GeV and tan/3 = 7. As we see, for large values of m 3 / 2 the relevant branching ratio 
exhibits a fixed-point behavior. The reason for this effect is that, while m\i% enters in the 
calculation via gaugino masses, m 3 / 2 multiplies the scalar matrices and thus dominates the 
flavor-violating processes in the case of non-zero off-diagonal contributions at the GUT 
scale. 

Figure |8| finally, shows the increase of the BR as tan /3 increases, for fixed values of 
m i/2 — 300 GeV and A = — 1.5m 3 / 2 (solid lines). Then, the experimental bounds impose 
severe constraints on the allowed range of (tan/3, myz, "7.3/2) that one may have. For 
certain combinations of these terms, we predict BR below the experimental limits. This 
means that the model we analyze, with non-universal soft terms at the GUT scale, is 
consistent with the current experimental limits for low masses of the sleptons and high 
masses for gauginos, when low values of tan/3 are considered (dashed line). 

Similar considerations may be done for r — > /ry. This is presented in fig.[|. Here, 
the current experimental bounds are not as strict as in the previous case. Then, in the 
framework of the models that we are discussing, we find that, for a wide region of the SUSY 
parameter space, the predicted BR is well below the experimental bounds. However, the 
prediction exceeds experimental limits for larger values of tan/3. 

Finally, in fig. |1(] we compare our results with the choice of parameters made in ||. 
The choice of 6 = 1/2 increases the values of the mixing terms of the scalar matrices at 
the GUT scale as can be seen from (|8|), hence for this choice of parameters the results 
obtained for the BR are one order of magnitude bigger than in our case. For completeness, 
we note that in || , the choice of the Yukawa coefficients for the lepton matrix ^ were 
en = 4.0, C12 = c 2 i = 0.9, c 22 = 1.08, c 33 = 1.9. In our case, we have instead c\\ = 1.0, 
C12 = c 2 i = 0.4, c 22 = 2.2, c 33 = 1.0 and 6=1. 
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In retrospect, we can infer from the figures that in the case of non- universality in the 
scalar sector the allowed ranges of mi/2, tan/5 and m 3 / 2 are extremely limited. This is 
rather evident in particular from our figure 8. Our conclusions are rather generic for this 
class of models, as far as mixing is also predicted in both the fermion and s-fermion mass 
textures. This fact naturally raises the question whether the simple [/(l)-models are ca- 
pable of generating a completely realistic low energy theory. We think that these problems 
cannot find a solution in the present models. Again, as stressed also in the introduction, 
flavor violations indicate that there should be a kind of 'mechanism' in the scalar sector to 
suppress large mixing effects. In our opinion, string derived models are probably the only 
realistic ones which may offer new 'mechanisms' of additional suppression. For example, 
in addition to the gauge and [/(l)-family symmetries, the trilinear and non-renormalizable 
terms in string derived models have to respect additional symmetries arising from mod- 
ular invariance constraints. As a result, a large portion of (L^l) -invariant) mixing terms 
are eliminated by these string symmetries. Additional suppression of the mixing effects 



in the scalar sector may also arise by cyclic symmetries as those discussed in ref|18 



In the Figures, we have shown the constraints we can obtain in m^, tan j3 and 7B3/2 from 
ji — > e7 decays. Let us finally see what this implies for the range of magnitudes of physical 
masses. To do so, we give the low energy sparticle masses for some indicative values of 
the input parameters: For a fixed value of tan (3 = 7, the lightest neutralino varies from 
80 to 196 GeV as m 1 / 2 scales between 200 to 450 GeV. For the same inputs, the lightest 
chargino varies from 196 to 365 GeV. In both cases the dependence of the results on m 3 / 2 
is small. For a fixed value mi/2 = 300 GeV, the masses of the lightest neutralino and 
chargino change by at most 3 GeV, while the masses of the heavier charginos/neutralinos 
decrease by up to 50 GeV, as tan (5 increases from 3 to 14. 

What about the charged slepton masses? Here, we find that changes with tan/? are 
negligible for m^j ~ 300 GeV. For tan (3 = 1 and m.1/2 = 200 GeV, when 7713/2 varies from 
100 to 500 GeV, we have: m fR » (125 - 500) GeV, m e - L » (160 - 520) GeV, while in this 
case the lightest scalar particle for our initial conditions is m^ ~ (100 — 500) GeV. For 
mi/2 = 450 GeV and the same range of m 3 / 2 , m? R « (200 — 520) GeV, mg L ~ (300 — 590) 
GeV, while the lightest scalar particle is again mi « (100 — 500) GeV. 



7 Conclusions 

In the present work, we investigated in detail the predictions of U(l) family symmetries 
for the rare processes fi — > ey and r — ► /ry. We worked in the small tan/? regime, and 
found that in this class of models, /x — > rare decays may occur at significant rates, 
particularly in the case of non-zero flavor mixing at the GUT scale. Demanding that the 
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predicted values do not exceed the current experimental limits, we can put important 
bounds on the model parameters. On the contrary, r — > /ry decays are less dangerous 
and thus do not lead to strong bounds. 

Focusing on /i — > ej, we first discussed the case where no mixing of soft terms at 
the GUT scale occurs. In this case lepton flavor violation arises from non-zero neutrino 
masses in the theory. We find that the dominant contributions to the decay rate cancel 
around a certain value of m.3/2, leading to very small values for the relevant branching 
ratio. An important observation is that the mixing effects of the Dirac matrix in the 
light sneutrino sector are canceled at first order. Here, we should stress that once the 
accuracy of the experiments is improved, it will be important to have a precise calculation 
which takes into account the precise form of the sneutrino mass matrix. The branching 
ratio of the process increases as m^j becomes smaller and tan/? larger, but we are still 
below the experimental bounds for typical values of the model parameters. Given that 
the experimental bounds will improve in the future, non-detection of fi — > e'j events will 
be associated to low values of gaugino masses as we can see in figure [|. 

The situation changes when mixing effects in the soft masses are introduced at the 
GUT scale. Since the effects of the off-diagonal terms are much larger than those of the 
massive neutrinos, the cancellation effects that we discussed are no longer present and 
large lepton number violating effects may be generated. As before, the branching ratio 
of the decay increases as becomes smaller and tan/3 larger, but now strong bounds 
on the model parameters are derived. In this scenario, the BR(fi — > ey) experimental 
bounds constrain the SUSY parameter space to large values for gaugino masses and small 
scalar masses as shown in figure §. 
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Figure 1: The generic Feynman diagrams for the \i — > cy decay. / stands for charged 
slepton (a) or sneutrino (b), while an d represent neutralinos and charginos re- 
spectively. 
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Figure 2: Amplitudes contributing to the BR(fi — > ey) when universal soft masses at 
the GUT scale are considered (A = 0). Dashed lines are the two contributions to A^. 
The curves are obtained using tan/3 = 7, m^j = 300 GeV, and A = — 1.5m 3 / 2 as input 
parameters. 
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Figure 3: BR(fi — > cy) for a range of values of mi/ 2 (labeled above). Universal soft masses 
at the GUT scale are considered (A = 0). The curves are obtained using tan/3 = 7 and 
A = — 1.5m 3 / 2 as input parameters. 
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Figure 4: BR(fi — > cy) for a range of values of tan/3 (labeled above). Universal soft masses 
at the GUT scale are considered (A = 0). Solid lines are obtained using m 1 / 2 = 300 GeV 
and A = — 1.5m 3 / 2 as input parameters. 
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Figure 5: BR(/j, — > cy) for a range of values of A (labeled above). Universal soft masses 
at the GUT scale are considered (A = 0). The curves are obtained using tan/3 = 7 and 
mi/2 = 300 GeV as input parameters. 




Figure 6: Amplitudes contributing to the BR(fi — > ey) when non universal soft masses 
at the GUT scale are considered (A ^ 0). The curves are obtained using tan/3 = 7, 
mi/2 = 300 GeV, and A = — 1.5m3/2 as input parameters. 
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Figure 7: BR(fi — > ej) for a range of values of m 1 / 2 (labeled above). Non universal soft 
masses at the GUT scale are considered (A ^ 0). The curves are obtained using tan (3 = 7 
and A = — 1.5m 3 / 2 as input parameters. 
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Figure 8: BR(fi — > ej) for a range of values of tan/3 (labeled above). Non universal 
soft masses at the GUT scale are considered (A ^ 0). Solid lines are obtained using, 
mi/2 = 300 GeV, and A = —l.hmz/2 as input parameters. 
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Figure 9: BR(t —>■ yrf). Non universal soft masses at the GUT scale are considered 
(A / 0). Experimental limits are violated for large tan j3 and small values of 777.3/2 (dashed 
line). In both cases A = — 1. 57773/2 ■ 
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Figure 10: BR(t — > fij), for = 300 GeV, tan/3 = 7, and A = — 1.5m 3 / 2 . The 
dashed line is obtained using the charge and coefficient choices of ||, while the solid line 
is derived using the choice of parameters presented in this paper. Non universal soft 
masses at the GUT scale are considered (A ^ 0) . 
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